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Abstract 



The COMPASS Collaboration at CERN has measured the transverse spin azimuthal asymmetry of 
charged hadrons produced in semi-inclusive deep inelastic scattering using a 160 GeV /i + beam 
and a transversely polarised NH3 target. The Collins asymmetry of the proton was extracted in the 
Bjorken x range 0.003 < x < 0.7. These new measurements confirm with higher accuracy previous 
measurements from the COMPASS and HERMES collaborations, which exhibit a definite effect in 
the valence quark region. The asymmetries for negative and positive hadrons are similar in magnitude 
and opposite in sign. They are compatible with model calculations in which the u-quark transversity 
is opposite in sign and somewhat larger than the d-quark transversity distribution function. The 
asymmetry is extracted as a function of Bjorken x, the relative hadron energy z and the hadron 
transverse momentum p'j. The high statistics and quality of the data also allow for more detailed 
investigations of the dependence on the kinematic variables. These studies confirm the leading-twist 
nature of the Collins asymmetry. 
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Deep inelastic lepton-nucleon scattering (DIS) measurements as a tool to unveil the structure of the 
nucleon started in the late 60's at the Stanford Linear Accelerator (SLAC), when for the first time a 
high energy electron accelerator became available. A wealth of eN scattering data was collected, and 
eventually it became clear that scattering at large momentum transfer could be interpreted as elastic 
scattering off the nucleon constituents, the "partons". From the dependence of the cross-section on the 
energy and the momentum transfered to the nucleon it was possible to identify the charged partons with 
the earlier postulated quarks. In the subsequent years the SLAC energy was gradually increased. Higher 
energy experiments could be performed at CERN and FNAL using muon beams from % and K decays, 
and at the HERA ep collider. The high beam energies allowed measurements at larger values of Q 2 , 
the square of the four-momentum transfered to the nucleon, and at smaller values of x, the fraction of 
the nucleon momentum carried by the parton. All these data in combination with the neutrino-nucleon 
data eventually allowed the extraction of the nucleon parton distribution functions (PDFs) and their Q 2 
dependence, in particular of the functions q(x, Q 2 ) (or ff(x, Q 2 )), which describe the number density of 
quarks of flavour q. 

In a second generation of DIS experiments, polarised lepton beams and polarised targets were used. The 
goal of these experiments was the measurement of the helicity distributions Aq(x) (or g\ {%)), defined as 
difference between the number densities of quarks of flavor q with helicity equal or opposite to that of the 
nucleon, for a nucleon polarised along its direction of motion (longitudinal polarisation). Following the 
discovery by the EMC at CERN in 1988 that the quark spins contribute only little to the proton spin, the 
interest in the nucleon spin structure was revived. In the past decade, a third generation of experiments 
at CERN (COMPASS), DESY (HERMES) and JLab began to study the nucleon spin structure using 
semi-inclusive DIS (SIDIS) processes, in particular the gluon contribution to the nucleon spin. These 
experiments also use polarised beams and polarised targets, while additionally hadrons produced by 
fragmentation of the struck quark are reconstructed and identified. 

Since the 90's |0]] it is well known that in order to fully specify the quark structure of the nucleon at twist- 
two level in quantum chromodynamics (QCD), the transversity distributions Ajq{x) (or h\(x)) have to 
be introduced in addition to the momentum distributions q(x) and the helicity distributions Aq(x). For 
a given quark flavour q, Ajq{x) is the analog of the helicity distribution in the case of a transversely 
polarised nucleon. Helicity and transversity distributions coincide in the non-relativistic quark model 
but they are expected to be different when relativistic effects are taken into account. As there exists no 
gluon transversity distribution, the Q 2 evolution of Ajq(x) is quite different from that of Aq(x). The first 
moment of the valence quark transversity distribution is related to the tensor charge 8q, 



Jo 

which together with the vector and axial charge characterises the nucleon as a whole. The tensor charge 
is presently being calculated with steadily increasing accuracy by lattice QCD (2J]. 

The transversity PDF is chiral-odd and thus not directly observable in inclusive deep inelastic lepton- 
nucleon scattering. In 1993 Collins suggested that it could be measured in SIDIS processes, where it 
appears coupled with another chiral-odd function Ji]], which by now is known as "Collins fragmentation 
function" A^D q (or H^ q ). It is the transverse-spin dependent part of the standard fragmentation function 
(FF) D q that describes the correlation of quark (q) transverse polarisation and hadron (h) transverse 
momentum. This mechanism leads to a left-right asymmetry in the distribution of hadrons produced in 
the fragmentation of transversely polarised quarks, which in SIDIS shows up as an azimuthal transverse 
spin asymmetry AqoU (the "Collins asymmetry") in the distribution of produced hadrons. At leading 
order this asymmetry can be written as 




(1) 



Acoii 



^ q e 2 q -A T q®A»D h q 
Z q e 2 -q®D h q 



(2) 
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where <g) indicates the convolutions over transverse momenta. The Collins asymmetry is accessed 
through the amplitude of the sin<J>c modulation in the hadron azimuthal distribution. Here the Collins an- 
gle <t>c = <j>h + $s — n is the sum of the azimuthal angles of the hadron transverse momentum pj (0/,) and 
of the spin direction of the target nucleon (<j> s ) with respect to the lepton scattering plane, in a reference 
system in which the z axis is the virtual-photon direction. 

A non-zero Collins asymmetry for the proton was first observed in 2004 by HERMES fl] using an 
electron beam of energy 27.6 GeV. These results provided first evidence that both Collins fragmentation 
and transversity functions are non-vanishing, although room was left to possible explanations of the 
observed signal in terms of higher-twist effects. Independent evidence of a non-zero and sizeable Collins 
function came soon after from measurements by the Belle Collaboration of the correlation between 
the azimuthal angles of hadrons in two jets resulting from e + e~ annihilations into hadrons, recently 
confirmed also by the Babar Collaboration |7J]. For a comprehensive review of recent experiments and 
theoretical developments see e.g. refs. ii and & 

Using a 160 GeV \i + beam COMPASS measured SIDIS on a transversely polarised 6 LiD target in 2002, 
2003 and 2004. In those data no size-able Collins asymmetry was observed within the accuracy of 
the measurements 110j4l2j], a fact that was understood in terms of a cancellation between the u- and 
d-quark contributions 111 ill . The COMPASS data are still today the only SIDIS data ever taken on a 
transversely polarised deuteron target and provide constraints on the d-quark contribution Q. Together 



with the HERMES and Belle data they allowed for the first global analyses and first extractions II 141. II 511 
of the transversity distributions for u- and d-quarks and of favoured and unfavoured Collins FFs. It is 
important to note that global analyses of e + e~ — > hadrons and SIDIS data are a necessity, since neither 
can the Collins FF be extracted from e + e~ data alone, nor can the Collins and transversity functions 
be disentangled using only SIDIS data. Similar considerations hold for the transverse spin asymmetry 
in hadron pair production in DIS lHoTl . where a signal on the proton was recently measured by HER- 



MES [17] and COMPASS [18]. In this case the transversity distribution is coupled with the di-hadron 



FF which was recently accessed lll9ll using the Belle data. 



In 2007 COMPASS measured for the first time SIDIS on a transversely polarised proton (NH3) target 
using the n + beam of energy 160 GeV, thus extending the measured x range to values about 10 times 
smaller as compared to HERMES. The results [20] for the Collins asymmetry in the valence quark region 
were in agreement with those of HERMES 12 ill , in spite of the considerably larger Q 2 values. 

As the Collins asymmetry includes convolutions of distribution and fragmentation functions depending 
on different variables, a large data sample is necessary in order to perform a multi-dimensional anal- 
ysis and test the properties of the observable, namely factorisation and evolution. At COMPASS, the 
detailed investigation of the kinematic dependence of the Collins asymmetry was a strong motivation 
to improve on the precision of the 2007 measurement, and hence the entire 2010 run was dedicated to 
SIDIS measurements using the transversely polarised proton target (NH3). 

In this Letter, first results for the Collins asymmetry from the 2010 data are presented. The Sivers 
asymmetry, which is the amplitude of a sin(<^ — (j) s ) modulation in the hadron azimuthal distribution, 



was extracted from the same data and the results are given in a parallel paper Il22ll . 



The COMPASS spectrometer J2J is in operation in the SPS North Area of CERN since 2002. The 
principle of the measurements and the data analysis were already described in refs. fills Eft!. In 



2010, the spectrometer configuration was very similar to that used in 2007 II20|] . Additionally, a new 
triggering system for large-angle muons was used, which is based on two large area scintillator counter 
hodoscopes with 32 horizontal bars each and a suitable coincidence matrix to provide target pointing in 
the non-bending vertical plane. The target is polarised along the vertical direction and consists of three 



The Collins asymmetry was very recently measured for the first time on a transversely polarised 3 He target at JLab flUl . 
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cylindrical cells with a diameter of 4 cm; the central cell is 60 cm long, and the two outer ones are 30 cm 
long and 5 cm apart. Neighbouring cells are polarised in opposite directions, so that data for both spin 
directions are recorded at the same time. 

The jit + beam had a nominal momentum of 160 GeV/c with a momentum spread Ap/p = ±5% and a 
longitudinal polarisation of -80%. The data were taken at a mean beam intensity of 3.5 • 10 8 jU/spill, for 
a spill length of about 10 s every 40 s. About 37 • 10 9 events, corresponding to 1.9 PB of data, were 
collected in twelve separate periods. In order to minimise systematics, in each period after 4-5 days of 
data taking a polarisation reversal requiring 1.5 days was performed and then data taking continued for a 
corresponding number of days. 

In the data analysis, in order to ensure the DIS regime, only events with photon virtuality Q 2 > 1 
(GeV/c) 2 , fractional energy of the virtual photon 0.1 < y < 0.9, and mass of the hadronic final state 
system W > 5 GeV/c 2 are considered, leading to a total of 16 • 10 7 DIS events. The charged hadrons are 
required to have at least 0. 1 GeV/c transverse momentum p\ with respect to the virtual photon direction 
and a fraction of the available energy z > 0.2. After these cuts about 8 • 10 7 hadrons are left and used for 
the extraction of the asymmetries. This is referred to as "standard sample" in the following. 

The asymmetries are measured separately for positive and negative hadrons as function of x, z or p\. The 
{Q 2 } values are about 3.7 (GeV/c) 2 in all z and p h T bins, while they strongly vary with x due to the fixed 
target kinematics (see fig. CD)- The mean values of x in the z. and p\ bins, of z in the x and p\ bins, and of 
Pj in the x and z bins are (x) = 0.05, (z) = 0.38, and (p h T ) = 0.52 GeV/c respectively. 

In every bin of x, z or p\, for each period of data taking the asymmetries are extracted from the number 
of hadrons produced in each cell for the two directions of the target polarisation. Using an extended 
Unbinned Maximum Likelihood (UML) estimator all 8 azimuthal modulations expected in the 
transverse spin dependent part of the SIDIS cross-section Il24n are fitted simultaneously. 



The measured amplitude of the modulation in 
sin<I> c is e c = D NN fP T A C oiu where D NN = (1 - 
y)/(l —y + y 2 /2) is the transverse-spin-transfer 
coefficient from target quark to struck quark, / 
the dilution factor of the NH3 material, and Pj the 
magnitude of the proton polarisation. In order to 
extract Ac u, the measured amplitudes Sc in each 
period are divided by /, Pj and D^n- The dilu- 
tion factor of the ammonia target is calculated for 
semi-inclusive reactions Il25ll and is evaluated in 
each x bin; it increases with x from 0.14 to 0.17, 
and it is assumed constant in z and p\. The pro- 
ton target polarisation (~ 0.8) was measured indi- 
vidually for each cell and each period. The final 
asymmetries are obtained by averaging the results 
of the 12 periods, after having verified their statis- 
tical compatibility. 
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Fig. 1 : Mean values of Q 2 in the x bins for standard 
sample range 0.1 < y < 0.9 (closed circles, •) and 
for events with 0.05 < y < 0.10 (closed squares, ■). 



Extensive studies were performed in order to as- 
sess the systematic uncertainties of the measured 

asymmetries, and it was found that the largest contribution is due to residual acceptance variations within 
the data taking periods. In order to quantify these effects, various types of false asymmetries are calcu- 
lated from the final data sample assuming wrong sign polarisation for the target cells. Moreover, the 
physical asymmetries are extracted splitting the events according to the detection of the scattered muon 
in the spectrometer (top vs bottom, left vs right). The differences between these physical asymmetries 
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Fig. 2: Correlation coefficient between the Collins and Sivers asymmetries as a function of x, z and p\. 



and the false asymmetries are used to quantify the overall systematic point-to-point uncertainties, which 
are evaluated to be 0.5 times the statistical uncertainties. For further details on the evaluation of the 
systematic uncertainties of the 2010 data see ref. [26]. A systematic scale uncertainty of 3% arises from 
the measurement of the target polarisation. 

The correlations between the Collins asymmetries obtained after binning the data in x, in z or in p\ were 
also evaluated. The Collins asymmetries measured along different projections of the (x, z, p\) phase 
space are statistically correlated, because their mean values have to be the same. In the COMPASS case 
these correlation coefficients are all smaller than about 0.2, but non-negligible, so that they should be 
taken into account in any global fit. 

The correlations between the amplitudes of the 8 modulations allowed in the UML fit are found to 
be small and about the same for positive and negative hadrons. No correlations are expected a priori, 
since these azimuthal modulations are independent, but in the COMPASS experiment they arise mostly 
because the coverage in the horizontal plane of the muon triggering system is not complete. Figure |2] 
shows the correlation coefficients between the Collins and Sivers asymmetries as a function of x, z and 
Pj, as given by the UML estimator. As can be seen, they are always smaller than 0.15. This correlation 
regards only the statistical fluctuations of the observables and is relevant only in case of simultaneous fits 
to the various asymmetries. 

Figure [3] shows the Collins asymmetries measured as a function of x, z, and p\ for positive and negative 
hadrons. The error bars are statistical only. As can be seen in the figure, the Collins asymmetry has a 
strong x dependence. It is compatible with zero in the small x region accessible at COMPASS and in- 
creases up to 0.05 in the valence quark region (x > 0.1), confirming that the trans versity PDF is a valence 
object. The data exhibit a mirror symmetry with respect to the hadron charge, which is interpreted as 
due to the same size and opposite sign of the favoured and unfavoured Collins FFs II 1011 . The values are 
in agreement with our previous measurements 112011 . with an important gain in statistics: the statistical 
and the systematic uncertainties are reduced by about a factor of two. In the same figure, the data are 



compared with the predictions of ref. Il27ll which were obtained by fitting the HERMES proton data |2i 
the COMPASS deuteron data Ill2ll and the Belle e + e~ fragmentation data |6J]. The observed agreement 
supports the weak Q 2 dependence of the Collins asymmetry assumed in the calculation. 



The high statistics of the 2010 data also allows for a thorough investigation of the kinematic dependences 
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of the Collins asymmetry. As a starting point, the z and p\ dependences of the asymmetries are also 
evaluated for the region x > 0.032, where the signal is different from zero. In this region the mean 
value of Q 2 is equal to 5.75 (GeV/c) 2 . As expected the results, given as open points in fig. |4j show a 
larger magnitude of the asymmetry due to the higher mean value of x, and are compatible with a linear 
dependence on pt up to about 1 GeV/c. These results agree in magnitude and sign with HERMES 
measurements [21] in this x region that is common to both experiments, although the mean values of Q 2 
in each x bin are 3 to 4 times larger in COMPASS. 

The Collins asymmetries are also studied at low y, namely in the region 0.05 <y < 0.1. Given the strong 
correlation between y and the mass of the hadronic final state system W, these events have W values 
between 3 and 5 GeV/c 2 , smaller than those of the standard sample. Also, as can be seen in fig. [El for the 
x range 0.032-0.70 {Q 2 ) is smaller by about a factor of 3 with respect to the events with 0.1 < y < 0.9. 
The measured Collins asymmetries for this low-y sample are given in fig. |5l where they are compared 
with the results obtained by splitting the standard y range into two bins, 0.1 < y < 0.2 and 0.2 < y < 0.9. 
The asymmetries measured as a function of x are all compatible and only give some hint for a decrease at 
higher y for negative hadrons. This is a further evidence that, if there is a Q 2 dependence of the Collins 
asymmetry, it has to be weak, which suggests both transversity and Collins functions being leading- 
twist quantities. Theoretical calculations of the Q 2 evolution of the transverse momentum dependent 
functions [29] are ongoing but not yet available for the Collins asymmetry. 

The z dependence of the Collins asymmetry is further studied extending the range towards the target 
fragmentation region, namely by measuring the asymmetries for hadrons with 0. 1 < z < 0.2 as a function 
of x and p\. The results are well compatible with those obtained from the standard z > 0.2 sample. 



All the results given in this Letter are available on HEPDATA 13011 . In particular, the asymmetries for the 
standard sample as functions of x, z and p\ have also been combined with the already published results 
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Fig. 4: Collins asymmetry as a function of z and p'f for positive (top) and negative (bottom) hadrons. 
The open points (o , slightly shifted horizontally) are the values obtained in the range x > 0.032. The 
closed points (•) refer to the complete x range and are the same as in fig. [3] 



from the 2007 run [20] and are also available on HEPDATA. 

In summary, this Letter presents the results from the 2010 COMPASS data for the Collins asymmetry 
on a proton target measured as a function of x, z and p\ using the entire collected statistics. They are 
compatible with the 2007 COMPASS results whereby improving the precision by a factor of two in 
the statistical and systematic uncertainties. The correlations between the Collins asymmetries measured 
along x, z and p\ respectively, and the correlations between the Collins and Sivers asymmetries have 
been investigated and found to be small. The Collins asymmetry has also been measured extending 
the kinematic range to smaller z and y values. The results confirm that the Collins asymmetry is a 
real leading-twist object. The availability of the precise COMPASS data on the proton, as well as the 
COMPASS data on deuteron, the HERMES proton results and the e + e~ fragmentation data, makes it 
mandatory to perform a new global analysis to disentangle the Collins and transversity functions and to 
study their properties. 

We acknowledge the support of the CERN management and staff, as well as the skills and efforts of the 
technicians of the collaborating institutes. 
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